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Abstract 

Applying effective Lagrangian method and on-shell scheme, we analyze the electroweak correc- 
tions to anomalous dipole moments of lepton from some special two loop diagrams in which a closed 
heavy fermion loop is attached to the virtual gauge bosons or Higgs fields. As the masses of vir- 
tual fermions in inner loop are much heavier than the electroweak scale, we verify the final results 
satisfying the decoupling theorem explicitly if the interactions among Higgs and heavy fermions 
do not contain the nondecoupling couplings. At the decoupling limit, we also present the leading 
corrections to lepton anomalous dipole moments from those two loop diagrams in some popular 
extensions of the standard model, such as the fourth generation, supersymmetry, universal extra 
dimension, and the littlest Higgs with T-parity. 
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I. INTRODUCTION 



At both aspects of experiment and theory, the magnetic dipole moments (MDMs) 
of leptons draw great attention of physicists because of their obvious importance. The 
anomalous dipole moments of lepton not only provide a potential window to detect new 
physics beyond the standard model (SM), but also can be used for testing loop effect in 
electroweak theories. The current experimental world average of the muon MDM is[l] 

a ex 'P = 11 659 208 ± 6 x 1(T 10 . (1) 

Contributions to the MDM of muon are generally divided into three sectors: QED loops, 
hadronic contributions as well as electroweak corrections. The largest uncertainty of the 
SM prediction originates from the evaluation of hadronic vacuum polarization and light- 
by-light corrections. Depending on which evaluation of hadronic vacuum polarization is 
chosen, the differences between the SM predictions and experimental result lie in the range 

a. 

For the electroweak corrections, the standard one loop contribution amounts to 19.5 x 
10~ 10 , and the one loop corrections from new physics sector are generally suppressed by 
A^ w /A^ p . Here A EW denotes the electroweak energy scale, and A NP denotes the energy scale 
of new physics. Comparing with the analysis at one loop level, the two loop analysis is 
more complicated and less advanced. Utilizing the heavy mass expansion approximation 
(HME) together with the projection operator method, the authors of Ref.j^] have evaluated 
the two loop standard electroweak corrections to muon MDM. Within the framework of CP 
conservation, Ref. [s| presents the supersymmetric corrections from some special two-loop 
diagrams where a closed chargino (neutralino) or scalar fermion loop is inserted into those 
two-Higgs-doublet one-loop diagrams. Ref.[6j discusses the contributions to muon MDM 
from the effective vertices H ± W Tr y ) /io(-^o)77 which are induced by the scalar quarks of the 
third generation. Furthermore, the contributions from two loop Bar-Zee-type diagrams to 
the electric dipole moments (EDMs) of light fermions are discussed extensively in literature 

In this paper, we calculate the corrections to the anomalous dipole moments of lepton 
from some special diagrams in which a closed heavy fermion loop is attached to the virtual 



electroweak gauge or Higgs fields. The effective Lagrangian method can yield the one loop 
electroweak corrections to lepton MDMs and EDMs exactly in the SM and beyond, and has 
been adopted to calculate the two loop supersymmetric corrections for the branching ratio of 
neutron EDM [9] and lepton MDMs and EDMs 10, 11]. In concrete calculation, 



S7 



we assume that all external leptons and photon are off-shell, then expand the amplitude of 
corresponding triangle diagrams according to the external momenta of leptons and photon. 
Using loop momentum translating invariance, we formulate the sum of amplitude from those 
triangle diagrams corresponding to same self energy in the form which explicitly satisfies 
the Ward identity required by the QED gauge symmetry, then get all dimension 6 operators 
together with their coefficients. After the equations of motion are applied to external leptons, 
higher dimensional operators, such as dimension 8 operators, also contribute to the lepton 
MDMs and EDMs in principle. However, the contributions of dimension 8 operators contain 
the additional suppression factor m|/A? comparing with that of dimension 6 operators, 
where mi is the mass of lepton. Setting A EW ~ lOOGeV, one obtains easily that this 
suppression factor is about 10~ 6 for the muon lepton. Under current experimental precision, 
it implies that the contributions of all higher dimension operators (D > 8) can be neglected 
safely. 

We adopt the naive dimensional regularization with the anticommuting 7 5 scheme, where 
there is no distinction between the first 4 dimensions and the remaining D — 4 dimensions. 
Since the bare effective Lagrangian contains the ultraviolet divergence which is induced 
by divergent subdiagrams, we give the renormalized results in the on-mass-shell scheme 



12j. Additional, we adopt the nonlinear gauge with £ = 1 for simplification 13J. This 



special gauge-fixing term guarantees explicit electromagnetic gauge invariance throughout 
the calculation, not just at the end because the choice of gauge-fixing term eliminates the 
, yW ± G Zf vertex in the Lagrangian. 

This paper is composed of the sections as follows. In section [Til we introduce the effective 
Lagrangian method and our notations. We will demonstrate how to obtain the identities 
among two loop integrals from the loop momentum translating invariance through an ex- 
ample, then obtain the corrections from the relevant diagrams to the lepton MDMs and 
EDMs. Section IIHI is devoted to the analysis and discussion in some concrete electroweak 
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models. In section [TV] we give our conclusion. Some tedious formulae are collected in the 
appendices. 

II. THE CORRECTIONS FROM THE RELATING DIAGRAMS 

The lepton MDMs and EDMs can actually be written as the operators 

jC = 6 a Ict^I F 

^ edm = -\ d t 1^1,1 . (2) 

Here, a^ u = i[7 M , 7^/2, / denotes the lepton fermion which is on-shell, F is the electromag- 
netic field strength, m [ is the lepton mass and e represents the electric charge, respectively. 

It is convenient to get the corrections from loop diagrams to lepton MDMs and EDMs 
in terms of the effective Lagrangian method, if the loop diagrams contain the virtual fields 
which are much heavier than the external lepton, i.e. m v 3> m l with m v denoting the mass 
scale of virtual fields. Since jj = j> = m l <C m v for on-shell leptons and jt — > m v 
for photon, we can expand the amplitude of corresponding triangle diagrams according to 
the external momenta of leptons and photon. The two loop diagrams also contain some 
virtual light freedoms generally, such as virtual neutrinos, charged leptons or photon, and 
it is unsuitable to expand the propagators of light freedoms in powers of external momenta 
obviously. In order to obtain the corrections to lepton MDM and EDM properly, we should 
firstly match the amplitude of two loop diagrams from full theory to that of corresponding 
diagrams from effective theory which is composed by the QED Lagrangian and some higher 
dimension operators of light freedoms, then extract the Wilson coefficients of those high 
dimension operators which are only depend on the masses of virtual heavy freedoms as well 
as the possible matching scales. Finally, we strictly analyze the amplitude of corresponding 
diagrams from effective theory to obtain the contributions from the virtual light freedoms 
to lepton MDMs and EDMs. As discussed in the section HJ it is enough to retain only those 
dimension 6 operators in later calculations: 
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Of = (Pj-fF ■ o^l 



°7 = {§^ lF -°^(wj) 



777 

eQ.m, - 

(3) 

with £> M = <9 M + ieA^ and c<j t = (1 =f 7s)/2. 

Certainly, all dimension 6 operators in Eq.(j3]) induce the effective couplings among pho- 
tons and leptons. The effective vertices with one external photon are written as 



Of = ^f 2 {((P + kf + P% + + t)l P j>}u 
ieO 



Of = ^%{(^ + ^)7p + 7p^Ho; 



or=^^7>,^- (4) 

If the full theory is invariant under the combined transformation of charge conjugation, 
parity and time reversal (CPT), the induced effective theory preserves the symmetry after 
the heavy freedoms are integrated out. The fact implies the Wilson coefficients of the 
operators Of 3 satisfying the relations 



where Cf (i = 1, 2, • • • , 6) represent the Wilson coefficients of the corresponding operators 
O t in the effective Lagrangian. After applying the equations of motion to external leptons, 



5 



we find that the concerned terms in the effective Lagrangian are transformed into 



=► (C+ + C 2 * + C+)Oi + (C+* + c 2 - + c+*)o 6 - 

= -^{^(C 2 + + C 2 ~* + C 6 + ) Jo"" J + i9(C£ + + C+) Z^Tfe j}i^ . (6) 



Here, 3?(- • •) denotes the operation to take the real part of a complex number, and 3(- • •) 
denotes the operation to take the imaginary part of a complex number. Applying Eq.(j2]) 
and Eq.([6]), we finally get 



In other words, the MDM of lepton is proportional to real part of the effective coupling 
C 2 + C 2 * + Cq , as well as the EDM of lepton is proportional to imaginary part of the 
effective coupling C 2 + C 2 * + Cq. 

After expanding the amplitude of corresponding triangle diagrams, we extract the Wilson 
coefficients of operators in Eq.Q which are formulated in the linear combinations of one 
and two loop vacuum integrals in momentum space, then obtain the MDMs and EDMs of 
leptons. Taking those diagrams in which a closed heavy fermion loop is inserted into the 
propagator of charged gauge boson as an example, we show in detail how to obtain the 
MDMs and EDMs of leptons through the effective Lagrangian method. 

A. The corrections from the diagrams where a closed heavy fermion loop is in- 
serted into the self energy of gauge boson 

In order to get the amplitude of the diagrams in FigHJ one can write the renormalizable 
interaction among the charged electroweak gauge boson and the heavy fermions F a ^ in 
a more universal form as 




(4tt) 2 



K(C 2 + + C7* + C 6 + ), 



L^c+ + C 2 * + C+). 



(7) 



WFF 



—W-'» 



(8) 



s 
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(b) 




FIG. 1: The relating two-loop diagrams in which a closed heavy fermion loop is attached to virtual 
bosons, where the diagrams (e,f,g) contribute the counter terms to cancel the ultraviolet 
divergence arisen by divergent subdiagrams in (a,b,c,d) respectively. 



where the concrete expressions of depend on the models employed in our calculation. 
The conservation of electric charge requires Qp — Q a = 1, where Q a ^ denote the electric 
charges of the heavy fermions F a ^ respectively. 

Applying Eq.([H]), we write firstly the amplitude of those two loop diagrams in FigfTJ For 
example, the amplitude of FigJT](a) can be formulated as 



x 



(qx-p- k) 2 - m 2 



{ie - g^aiZp + k- 2q 1 ) p + 2(g pil k a - g^k^) } 
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X 



(?i - P) 2 ~ m l (<?i - P) 2 ~ m l 



Tr 



Hg2-gi + ra Fa , 
(92 - 9i) 2 - ™ 2 



(^)>{c>-+0 



^ (92 - P) 2 - m 2 



(9) 



Here A RE denotes the renormalization scale that can take any value in the range from the 
electroweak scale A EW to the new physics scale A NP naturally, and we adopt the shortcut 
notations: c w = cos# w , s w = sin# w , with 8 w denoting the Weinberg angle. Additionally, 
p, k are the incoming momenta of lepton and photon fields, p is the Lorentz index of photon. 
Certainly, the amplitude does not depend on how to mark the momenta of virtual fields for 
the invariance of loop momentum translation. It can be checked easily that the sum of 
amplitude for diagrams in Fig{T] satisfies the Ward identity required by the QED gauge 
invariance 



k p A(p, k) = e[S ww (p + k) - S ww (p)] , 



(10) 



where -4 ww p denotes the sum of amplitudes for the diagrams (a), (b), (c) and (d) in FigO], 
as well as S ww denotes the amplitude of corresponding self energy diagram, respectively. 

According the external momenta of leptons and photon, we expand the amplitude in 
Eq.® as 



d D qi d D q 2 



1 



^ww a >' k ) V ' A - J (2n) D (2tt) d ql{ql - mlf{{q 2 - q,) 2 - m 2 )(q 2 - m\ ) 



f 2 9l -(3p + fc) , 2q x -p 2p 2 + (p + k) 

X < 1 H 5 5 I o 5 " - 

[ qf — mr q 2 — m 



P 



q 2 — m 2 



92 -m 2 



4fa ■ (p + k)f + 8(gi ■ p)(g l ■ (p + k)) + l%x ■ p) 2 A(q 2 -pf 



(qf -ml) 2 
4:( qi -(3p + k))(q 2 -p) 



(9! - m 



2 ^2 



{ql - m 2 ){ql - ml ) 



■0/ - 9^(2p + k - 2q 



Up 



+2((7p At A; - 9pakp 




Tr 
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X7,{C>- + C> + }(^2-V+^) 



(11) 



since we only consider the corrections to lepton MDM and EDM from dimension 6 operators. 
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Because the denominators of all terms are invariant under the reversal q\ — ► —qi,q 2 — ► 
—q 2 , those terms in odd powers of loop momenta can be abandoned, and the terms in even 
powers of loop momenta can be simplified by 

d D q 1 d D q 2 qirfiuqipqivqiaqip, qxvaxvqipqioqiafhp 



(2tt)d (2tt) d ((q 2 



qi) 2 - ™%)($ - m 2 )(q% - m 2 ) 



S 



pvpaafi 



d D qi d u q 



D, 



D(D + 2)(D + 4) J (2ir) D (2it) d ((q 2 
d D qi d D q 2 



qi) 2 - ml){ql - m 2 )(q% - m 2 ) 



9l/i9ly9lp9le7-92a92/3 



(2n) D (2n) D {(q 2 - Ql ) 2 - m 2 )(q 2 - m\){q 2 - m, 
d D qi d D q 2 1 
(2tt)^ (2n)° ((q 2 - q,) 2 - m^q 2 - m 2 )(q 2 ~ < 



x 



Dq 2 M, ■ q.Y - (qlfql 



-s 



q 2 Mi ■ q 2 



D(D - 1)(D + 2)(D + 4) 
d D qi d D q 2 



D(D-l)(D + 2) 

qip,qivqi P q2aq2pq28 



pupcr J ap 



x 



(2n) D (2n) D {(q 2 - Ql ) 2 - m*)(q 2 - m\){q 2 - m. 
d D qi d D q 2 1 
(2ir) D (2n)° {{q 2 - q,) 2 - m 2 )(q 2 - m 2 )(q 2 - < 

-(D + l)q 2 qi ■ q 2 q 2 2 - 2(q, ■ q 2 f (q, ■ q 2 f 

& as ~r 



q% ■ q 2 q 2 



D(D — 1)(D + 2)(D + 4) ^ pap5 D(D-l)(D + 2) 

9,s9 Pl3 ) + g^(g„ a g P s + g vS g po ) + g^(9 m 9 pP + g„ g Pa ] 



pa ^.9 vpQ p$ 



(12) 



and those similar formulae presented in Eq.(5) of Ref[8|], where the tensors are defined as 



r ^pvpa QpvQpa QppQva QpaQvp 5 

S ,= q T _ + q T . + q T a + q T a + q a T 

pvpaafl J ph> paap J pp voctfl J pa vpa.fi J pa vpofi ^ pp vpaa 



(13) 



Summing over those indices which appear both as superscripts and subscripts, we derive 
all possible dimension 6 operators in the momentum space together with their coefficients 
which are expressed in the linear combinations of one and two loop vacuum integrals. In 
a similar way, one obtains the amplitude of other diagrams. Before integrating with the 
loop momenta, we apply the loop momentum translating invariance to formulate the sum of 
those amplitude in explicitly QED gauge invariant form, then extract the Wilson coefficients 
of those dimension 6 operators listed in Eq.([3]). Actually, we can easily verify the equation 

d D Ql d D q 2 



qip. 

(2tt)^ (2tt)^ (q 2 - m 2 )(q 2 - m 2 ){{q 2 - Ql ) 2 - m 2 ) 



= 



(14) 
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Performing an infinitesimal translation q\ — > q\, q2 — > q2 — cl with a p 
one can write the left-handed side of above equation as 

d D qi d D q 2 



0(p = 0, D), 



II 



(2tc)° (2tt) d {q\ - m 2 )(g 2 - m 2 )((g 2 - q^ - rng) 
d D qi d D q 2 



(2n) D (2tt) d (qf - m\){ql - m|)((g 2 - gi ) 2 - rng) 



x 



2g 2 • a 2(q 2 - qi) ■ a 



+ 



q\ - m 2 (g 2 - 9i) 2 - ml 



+ 



}■ 



(15) 



This result implies 



// 



d D gi d D q 2 



9i • 92 



(2tt) d (2tt) d (gf - m?)(g 2 - m|) 2 ((g 2 - gi ) 2 - m 2 



<Fgi d u q 2 



9i - 9i ■ 92 



(16) 



(2tt)^ (2tt)^ (g 2 - m 2 )(g 2 - m 2 )((g 2 - gi ) 2 - m 2 ) 2 " 

In a similar way, other identities listed in Ref.|8] can be derived. Using the expression of 
two loop vacuum integral [14j 



A 



d D qi d u q 2 



(2^ (2tt)^ (g 2 - m 2 )(g 2 - m 2 )((g 2 - g x ) 2 - m 2 ) 

a 2 r»(i + c )^_ 1, 

\ — 2 V ' 1 ' 



•?-2 



2(4vr) 4 (1 

H — (^2(xo In xq + xi In x\ + x 2 In x 2 ) — xo — xi — x 2 
— 2(x + X\ + x 2 ) + 2(x lnx + X\ lnxi + x 2 lnx 2 ) 
— x In 2 x — X\ In 2 xi — x 2 In 2 x 2 — $(x , Xi, x 2 



(17) 



and 



$(x, y, 2) = (x + y — z) lnx In?/ + (x — y + z) lnxlnz 

+ (— x + y + z) In y In 2 + sign(A 2 )\/|A 2 |\l/(x, y, z) 



^— (x, ?/, 2) = lnx In 7/ + In x In z — In ylnz + 2 lnx + — — ^ ^ ^(x, y, 2) 
ax ' ' . / 1 A 2 1 



(18) 



one obtains easily 



A 4e a 

RE 



a 2 ax 



9? 



(2tt)* (2tt)^ (g 2 - m\){q\ - m 2 )((g 2 - g x ) 2 - rng) 
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A t r 9 d , r r r d D qi d D q 2 qi ■ q 2 



A 2 \dx dx 2 nJ J 



(2tt)^ (2tt)^ (q 2 - mlM - ml)((q 2 - Ql ) 2 - ml) 

= 2( ^ )4 ^l H 2 l - ^^4^ + 7M 1 + 21nxo) + 2x2(1 + lna;o + lnX2; 

— (xi + x 2 ) In 2 x — (xi + 2x 2 ) lnx In x 2 — x 2 In 2 x 2 — Xi lnx lnxi + X\ lnxi lnx 2 
-2(x 1 +x 2 )lnx - 2x 2 lnx 2 - Xl ^ X ° — i^i]/^ Xl X2 ) j (19) 

which is equivalent to the identity Eq.( jT6i) . Here, e = 2—D/2 with D denoting the dimension 
of space-time, A is a energy scale to define Xj = m 2 /A 2 , and x R = A 2 E /A 2 . Additionally, 
A 2 = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, and the concrete expression of \I/(x, y,z) can be 
found in the appendix. Actually, the equation Eq. ffT9l provides a crosscheck of Eq. ffTT|) and 
Eq. (fl8l) rather than a verification of Eq. (fl6|) . In the limit z <C x, y, we can expand $(x, 2) 
according z as 

z 2 2 3 
$(x, y, z) = (p Q (x, y) + ztpx(x, y) + -^^(x, y) + -^^(x, v) 

+2z(\nz - l)'K 1 (x J y) + 2z 2 (^- - j)7t 2 (x,y) 



2! 4' 

-2z »(—--) Vi {x,y) + ■■■ (20) 



with 



^i(x,y) = 1 + Q 1A {x,y), 

x + y 2xy y 

7i 2 (x,y) = — - Ti ln - ; 

[x — yy (x — yy x 

l \ ^ 12xy 6xy(x + y) y 

n ^ y) - ( x -y)2 ( x _j,)4 ( x _ y )5 111 x > 

and the concrete expressions of function </?j(x, y) (z = 0, 1, 2, 3) can be found in appendix. 

After applying those identities derived from loop momentum translating invariance, we 
formulate the sum of the amplitude for Figfjja,b,c,d) satisfying QED gauge invariance and 
CPT symmetry explicitly, and extract the Wilson coefficients of those operators in Eq.Q. 

Since only the operators Of 3 actually contribute to the MDMs and EDMs of leptons 
when the equations of motion are applied to the incoming and outgoing leptons separately, 
the relevant terms in the effective Lagrangian are formulated as: 

ff _ (4tt)V rd^d^ 1 

stQ f RE J W D W D qM - ™ 2 ) 2 (fe - <?i) 2 - m 2 F J(q 2 - m 2 ^) 
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f^L^R , >fl*> L W(l) , />L*> L _ £R*£R W(2) 

+m m ± (c L ;c* + c R ;c L J^ (3) fe~ + a- 

X \ a P a P a P a/3/ ww V 2 3 



+ •••, (22) 

where Q f — — 1 represents the charge of leptons, and the expressions of form factors Af® (i = 
1, 2, 3, 4) are presented in appendix. 

Integrating over loop momenta, one gets the following terms in the effective Lagrangian: 

.T 2 (l+e) ' 



£eff = V2G F a e x w 



ns 2 w Q f {l-e) 2 



X 



5 x f„ + X I 



24e 



+ 24x 2 ^'^'^V 4 



x z 



or 



5 , 
H mx 



27 24 



9x„ 



+ r 1 (x w ,x_ ,x_ ) (e>r + er) 



£l,l (^Fc 5 -^Fo ) 



+ (C L *C R +C R *C L )(x x W 2 T— 1 

^V S «/3 S 'a/3 ^ S a/3 ^ cc0 J ^ F f} > Y\2sX 2 YlX 

w w 

19 5 i 
+ 72^ " 12^ ln X - + T *^ x *~ ' ^ + ^ 

w w 

+(CC - CC)(\^)^(^\.s)](^ - + 



(23) 



where G F = 1.16639 x 1CT 5 GeV~ 2 is the 4-fermion coupling, and a e = e 2 /4ir. Note that 
the above result does not depend on the concrete choice of energy scale A, and the concrete 
expressions of Ti(x,y,z), g id (x,y) (i, j = 1, 2 • • •) can be found in appendix. 

The charged gauge boson self energy composed of a closed heavy fermion loop induces the 
ultraviolet divergence in the Wilson coefficients of effective Lagrangian, the unrenormalized 
W ± self energy is generally expressed as 



(p, A RE ) = A 2 A- 9flu + (A? + % i AZ + ..-) (p V - PlJh ,) 



A 2 ' 

+ (B 1 W + ^5 2 W + ---)^, 



(24) 



where the form factors A^ l2 and Bf 2 only depend on the virtual field masses and renor- 
malization scale. Here, we omit those terms which are strongly suppressed at the limit of 
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heavy virtual fermion masses. The corresponding counter terms are given as 



Z^(p, A re ) = -[SmKA^) + m 2 JZ w (A nB )]g, v - 5Z w (A RE ) 
The renormalized self energy is given by 

£>, A RB ) = S-(p,A RE ) + S; c (p,A RE ) . 
For on-shell external gauge boson W , we have 



V 9nv ~ V^Vv 



(25) 



(26) 



12 



1 







lim 



E r,(^ m w)^(p) = e M (p) 



(27) 



w i w 

where e(p) is the polarization vector of W ± gauge boson. Inserting Eq. ff2^|) and Eq. fT25|) 
into Eq. (1271) . we derive the counter terms for the W ± self energy in on-shell scheme as 

,2 



m 



5Z os (m ) = AJ + -^-A? = A™ + x AI 

w Vw/ i A 2 i wz 

(5m 2 ' os (m ) = A" A 2 - m 2 5Z os 



(28) 



We should derive the counter term for the vertex jW + W~ here since the corresponding 
coupling is not zero at tree level. In the nonlinear gauge with £ = 1, the counter term 
for the vertex / yW + W~ is 



i8C lW + w - = ie ■ 5Z n 



( A re) [g^uih - k 2 ) p + g up (k 2 - h)^ + g Pi i{h ~ h)u 



(29) 



where hi (i = 1, 2, 3) denote the incoming momenta of and photon, and p, u, p denote 
the corresponding Lorentz indices respectively. 

We present the counter term diagrams to cancel the ultraviolet divergence contained by 
the bare effective Lagrangian in Fig{T|e,f,g), and we can verify that the sum of corresponding 
amplitude satisfies the Ward identity required by the QED gauge invariance obviously. 

Accordingly, the effective Lagrangian from the counter term diagrams is written as 

- 2 ; r(l + e) r ™r 5 We he 



5C 



c 



2s 2 A 2 Q 

he 
+ 12x" 



■ (Attx i 



12x 2 



72x 2 12x 2 



In a; 



{A- + x v A-)}(o; + o s 
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V2G 



5 



-(47TX n ) 



f 5 X F Q + ^Fo 



24e 



24x 2 

w w 

+ (C L *C R +C R *C L )(x x ) 1/2 



V27 + 24 



x 



9x., 



(o: + o; 



72x 2 + 12x 2 



\nx 



12ex 2 

w 

+ ••• 



12 X 2 ^ 1 ' 1 ^ Xp <=' ' Xp p' 



(30) 



Adding the counter terms to bare Lagrangian Eq. (123]) , we cancel the ultraviolet divergence 
there. However, the diagrams in Fig{TJ include the virtual neutrino which belongs to light 
freedoms contained by the effective Lagrangian. It is unreasonable obviously in the above 
analysis that the propagators of virtual neutrino are expanded according to the external 
momenta. In order to obtain the corrections to lepton MDMs and EDMs from light freedoms 
properly, we match the sum of amplitude from full theory to that from effective theory 15] 
at first: 



(31) 



where A MA represents the matching scale, and m Vh , m vl denote the masses of virtual heavy 
and light freedoms respectively. The left-handed side of above equation denotes the am- 
plitude from Figfj] derived through the above steps, the first term of right-handed side 
is the corrections to effective Lagrangian from heavy freedoms only, and the second term 
of right-handed side is the corrections obtained unsuitably to effective Lagrangian from 
light freedoms. Shrinking those heavy freedoms of Figfj] in a point, one obtains the cor- 




FIG. 2: The diagram of effective theory corresponds to those diagrams in Figfj] 



responding Feynman diagrams of effective theory in Figj2j Expanding the amplitude of 
diagram for effective theory in powers of external momenta, we can derive the corrections 
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J2 f^.{ m vn ^ma)^ T which are originated from light freedoms purely. Inserting the concrete 
expressions of f^,(m vl , A MA ) into Eq.( l3Ti) . one gets the corrections f^_(m Vh , A MA ) which are 
originated from heavy freedoms only Finally we analyze the amplitude of FigfS] strictly to 
get the corrections from virtual light freedoms to lepton MDMs and EDMs. Because the 
effective coupling among leptons and photon in Figj2] is induced by the dimension 8 opera- 
tors at least, the corrections from Figj2]to lepton MDMs and EDMs contain the additional 
suppression factor m^/A^ w comparing with that only from the heavy freedoms. Under our 
approximation, the resulted lepton MDMs and EDMs are respectively formulated as 

V w 

+2(x Fa x F0 ) 1/2 3* )T 3 (x w , x Fa , Xrp )}, 
d%F = -fj|^xJx FQ ^)V^(C^)T 4 (x w ,x FQ ,^) , (32) 

V w 

which only depend on the masses of virtual fields. It should be clarified that the correc- 
tions to lepton EDMs from the diagrams (a,b,c,d) in Figfj] do not depend on the concrete 
renormalization scheme adopted here since the relevant terms from bare Lagrangian do not 
contain the ultraviolet divergence. Using the expansion of <&(x,y,z) in Eq( f20l) . we get the 
asymptotic expressions of Ti(x, z,u), (i = 1, • • • , 4) at the limit z, «>xas 

rp ( \ 3 - 3Q/3 3 - 3 - 2Q/3 , N 

T^x.z.u) ~ + — Inu n^u) 

ox ox ox 

(3-2Q /3 )z + 6udir 1/ , (3 - 2Q p )zu + (3 + 2Q p )u 2 d 2 ^ . N 
8x ou 8x Ou 2 

1 2 <9 3 vr 2(4-Q )z-(l-2Q )u 

+ l2~x U (^-«)^.«) T7f x ^ M 

7 . ,dir 2 , . (z — u) 2 , . 

+ ife n( "- n) ^7 (z ' n) + ^4— ^ Z ^ + --- ' 

m/ , lnu 1 (5 — Qb)u dix, . . u(z — u)d 2 ix,. 

t 2 x, z, u) ~ -— + — tt, (z, u) - 1 , — -jriz, u) + -4 — -inr & «+••-, 

8x 8x Ax Ou 8x Ou 2 

1 - 3Q/3 1 + Q p dir l-Qpdu u{z- u) <9 3 tt 1 

T3(X ' M) " "W + ~^x~^ U) + -^x~^ U) ~ -^6x~^ U) 
Q,z - (3 + Q,)u _ 1-Qf, _ 8^ 



x du 2 ' 8x dzdu 

3 , , 1 — Qn , , rdn n dn n 



+ ^ M '^6X~ {Z - U) 
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' 2 



+ 



<9w cfe 



z.u) + 



T 4 (x, z,u)~-- — + u) - — — {z, u) 

8xu 8x ou 8x az 

This implies that the leading contributions contained in the asymptotic form of Eqj32] under 
the assumption m F = m Fa = ^> m w can be written as: 

° ; ™ K l^{ (18% - 13) ( |c i |2+| Ci 2 ) 

w 

+3(Q, - 3)(|^| 2 - lO 2 ) + ll»Oi)} + • " " , 
G F a e em ; (2 + Q /3 ) fl L 
^ 16V2^V + ' ( j 

v w 

where ellipses represent those relatively unimportant corrections. 

Comparing the result in Eq. fl32l) . the contributions from the corresponding diagrams 
contain the additional suppressed factor m 2 /A 2 5W when both of virtual charged gauge bosons 
in FigfT|a,b,c,d) are replaced with the charged Goldstone G^. However, we should consider 
the corrections from those two loop diagrams in which one of virtual charged gauge bosons 
is replaced with the charged Goldstone since it represents the longitudinal component 
of charged gauge boson in nonlinear gauge. For many extensions of the SM contain the 
charged Higgs, we also generalize the result directly to the diagrams in which a closed heavy 
loop is attached to the virtual if and fields simultaneously. 



B. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual W ± , G ± (F ± ) b osons 

Similarly, the renormalizable interaction among the electroweak charged Gold- 
stone/Higgs G ± (H^ and the heavy fermions F a ^ can be expressed in a more universal 
form as 



C 



S±FF c 

<->x* 



G-F a (g c a fu„ + ggu+)F p + H-F a (H^u. + H c >«u + )Fp + h.c. , (35) 



where the concrete expressions of G^f ,R , depend on the models employed in our 

calculation, the conservation of electric charge requires Qp—Q a = 1- Generally, the couplings 
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(a) 




FIG. 3: The two- loop diagrams in which a closed heavy fermion loop is attached to the virtual 
W ± , G ± or #± bosons. In concrete calculation, the contributions from those mirror diagrams 
should be included also. 



among the charged Goldstone/Higgs and leptons are written as 

em, 



V2> 



m s 



G l<jj-V l + B C H lu)-V l + h.c. 



(36) 



where the parameter B c also depends on the concrete models adopted in our analysis. In full 
theory, the couplings in Eq. fl35l) induce the corrections to lepton MDMs and EDMs through 
the diagrams in FigJ3J and the corresponding diagram of effective theory is same as that 
presented in FigJ2J 

After the steps taken in WW sector, the corresponding corrections from the diagrams in 
Figj3] to the lepton MDMs and EDMs are formulated as 



wg _ Gr F a e m [ 



32-\/27r-.s 



x„ 



1/2 



+ (^) / f 2 (x ,x ,x F ,x F m(g c ' L ( R +g c ' R c L 

1 V n- I z Vw) wJ F a 1 Fg) \ 0a ^ 0a ^>af 

+ (^y/ 2 F j x )$t(g c > L ( R -g c > R ( L )\ 

\ T / 4 V w> w) F a F (j ) V /3a a/3 Pa ' a0 J I 



d wo = G F a e em l ^ 



64^2 



3 e 2 



TC°S 
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,WH 



+ (^) 1/2 ^(x w ,a; w ,a; Fa ,x F/3 )SfeX L , " C 

+(^) 1/ V 4( ,.,,.,^,x FB )9fe'C - e;fci)} . 

w ^ w 

+ (^k) 1/ V 2 (a; w ,x ± ,x F ,x F )3?fH c ; i C / l+^ c ; R C L ,) 

\™ J -^V W H ± ' Fq, ' Fp) \ pa^ap pa ^ a/3 J 

\ j. / -5V W 5 H ± 5 F a 1 F^/ ^ /3q ^ Q/ 3 pa^apj 

\ I 2 



F 4 (x w ,x ,,x F ,x F ^(V/C* - K' R ( L „) \ , 

*\ w! H ± ' F a ' Fq) \ pa^ap /3a ^a/3J ' 



U>1 — 7= „ " X„ 



3„2 



7l°S 



+ f^k) 1/2 F 2 (a; ,x + ,x F ,x F )%(h c ' l ( r + H c ' r ( l ) 

\ «. / ZV w ' H± ' F a' */3' \ pa^ap pa^apj 

)sf^ c ' L c L - n c ' R c R 



/!„ \ 1/2 



,->-» ,->-» rp 

X W ± ' F a ' Fq 



)%(h c ' l ( r -n c ' R c L ) 

I \ pa ^ap pa ^ap J 



(37) 



The expressions of form factors y, 2, (i — 1, • • • , 4) can be found in appendix. 

Using the asymptotic formulae of form factors F iy (i = 1, •••,4) under the condition 

z, u > x, y, 



Fi{x,y,z,u) 



3(2 - 3Qp)z 3 - 32z 2 u - (20 - 9Qp)zu 2 , : 



3(z -uY 

(11 - 54Qp)z 2 - (151 - 54Qfi)zu + 2u 2 
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r d d 2 



9(z-u) 3 
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u 



2Q hl (x,y)n l (z,u) + ip 1 (z,u) 

d_ 

du 

d 2 1 r 1 

-m(^ J |0 M (a;,j/)7r 2 (2,u) + -ip 2 {z,u) 



+ 



F 2 (x,y,z,u) 



3z 3 - (35 - 2Q p )z i u - (29 - \2Q p )zu 2 + 15(1 - Q p )u 3 z 



3(z-uY 



u 
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F 3 (x,y,z,u) 



F±(x,y,z,u) 



(11 - 18Q p )z 2 - (223 - 90Qp)zu + (74 - 72Q p )u 2 
9(z -u) 3 
d d d 2 iv 

— + (1 - Qp) — - u—\ [2g 1A (x, y)ir 1 (z, u) + (fi(z, u) 

+ [(2 - Q p ) + (z + 9u)— + (1 - Qp)(z - u) — 

d 2 i r 1 
-u{z - u)—\ [e M (a;,y)7r a (z,u) + -<P2{z,u) 



10(z-u) r 



+ 



Qi A (x,y)7c 3 (z,u) + -y> 3 (z, w) 

(l + 4Qp)z 2 + (5-4Qp)zu z 
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(z — u) 2 

d r 

+ (1 + 2Q f3 ) — \2g 11 (x,y)n 1 (z,u) + <fi(z,u) 
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9 <9 n r 

^ + (1 [ 2 £i,i( x >^) 7r i(^ M ) + <Pi(z> u ) 

d 8 I r 

Qp-(z- u)— - (1 - Q^z - ^)t^J ^ ia (x,?/)7r 2 (2;,u) 



+ -<^ 2 (z,u) 



+ ■ 



simplify the expressions of Eg. (I37j) in the limit m 



f m Fa = m F« > m w as: 



G g g e mfm w 
32 v / 27r 3 s 2 m p 



9 _ Hg + (3 + 9l) ln I^l ^^ Cl z C L + gc,R C R \ 

4 18 3 m 2 -I V /3« S Q/3 /?« S "/3/ 
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2(9 -AQ ) 2(3 -Qfs) ln m 2 F 
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WG 
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2(5 - 9Qp) _ 2jl + 3gg) ^ m?_ 
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- ^ + (3 + %(lnm^-^K,m^) 



\ pa J a:p pa 3 ap 



+ 



+ 



19 - 20Q^ 2 - 4Q /3 



9 3 
16 _ 2 + 6Q^ 

9 3 



lnm 2 - ftl (m 2 ,m2 ))l3fjf^ c ' i C fl + H C ' R ( L 



(in 



V /3a 'a/3 pa 'a/3 



2^_6-2Q^ 2 (m 2 2 j^/^y _ H c,r c 



G F o; e em ; m w i3 c f r21 5 



ii 3 s 2 m. 



-^ + (3 + |)(K- ft!l ( ffl >^)) 



64^ 

\ /3a 'a/3 /3a '< 

19 - 20Q^ 2 - 4Q /3 



x 
+ 
+ 
+ 



9 3 
16 2 + 
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(in 



m 



w' £f±Vj V /9« afl /5a 'a/3 



(in 



^MV w' H±' J\ \ Pa ^ap pa ', 



2Q^_ 6 2Q^ / 2 _ 2 2 } xi / C , L * _^ c lO (3g) 

Q g V F 1,1 V w H±'/\ \ Pa^aP pa^apji V / 



The results indicate that the corrections to a/, di from the diagrams in Fig J3] are suppressed 
m F = m F 3> m w unless the couplings ?{ C > L > R violate the decoupling 



in the limit m I 



theorem. 



C. The corrections from the diagram where a closed heavy fermion loop is inserted 
into the self energy of Z gauge boson 

In order to get the amplitude of the diagram in FigJU we write the renormalizable 
interaction among the electroweak neutral gauge boson Z and the heavy fermions F a fi in a 
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(c) (d) 



FIG. 4: The two loop diagrams in which a closed heavy fermion loop is inserted into the propagator 
of virtual neutral gauge boson, where diagram (b) contributes the counter terms to cancel the 
ultraviolet divergence arisen by divergent subdiagram in (a). 



more universal form as 

£ ZFF = ^V ZM ^(C>- + , (40) 

w w 

where the expressions of depend on the concrete models employed in our calculation, 
and the CPT symmetry requires £^ = ^ = 




(■) (b) 



FIG. 5: The diagrams of effective theory correspond to those diagrams in FigUJ Where the diagram 
(a) corresponds to FigHJ^a) and (b) in full theory, and the diagram (b) corresponds to FigSJc) and 
(d) in full theory respectively. 



It is easy to check that the divergent amplitude of FigJU^a) satisfies the Ward identity: 

k p A$ p (p, k) = e[E_(p + k) - KM] , (41) 
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where p, k are the incoming momenta of lepton and photon fields, p is the Lorentz index of 
photon, A^' p denotes the sum of amplitude for the diagrams in Figj5](a), and S zz denotes the 
amplitude of corresponding self energy diagram, respectively. After a tedious calculation, 
the bare Lagrangian from the diagram FigJD^a) can be written as 

' 1 X F r , ' X Fa 



£eff = 



x- 



X F a + X Fp /4 



X, 



V3 



Q2,l( X F a , X F f3 ) 2 , . 

x 2 3x ' Xpa ' Xpfi 



x 



{Tf - Qf s i) 2 (°2 + o 8 -) + Qyjo: + o 3 - 
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K X F a X Bg. 



a/2 



2 (x Fa , x F ) 



lnx 



£X Z 



X z 



3x 2 



r ;-0/O (o a -+o.-)+o;<(o : 



^^(^Fo'^F.) X F a X F g .3 

^ ^ — - — \-inx, 



x c 



1 X Frv X I 



-QQ f U L «C +^)s 2 (Tf -Q f s 2 )( x F x f ) 1/2 

@i,i(. x F a i x f«) 1 + lnx, 



x z 



x z 



(42) 



where T z = — 1/2 represents the isospin of charged leptons. 



Generally, the unrenormalized self energy of the weak gauge boson Z can be written as 



fa A re) = A 2 A^ + (A? + + ■ • • ) (p V - PmPi/ ) 



A 2 ' 

+ (Bl + ^ 2 B z 2 + ...)p,p u . 



(43) 



Here, we omit those terms which are strongly suppressed at the limit of the large virtual 
fermion masses, and those form factors are actually decided by the masses of virtual fields 
and the renormalization scale. Correspondingly, the counter terms are given as 



= " ^!(A RE )+m 2 5Z(A RI 
The renormalized self energy is expressed as 

^ ( P ,A aB ) = E^( P ,A aB ) + Ej(p 1 A 
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V 9»v ~ V^Vv 



(44) 



(45) 



For on-shell external gauge boson Z, we have 12] 



1 



lim 



K >,mje» = e » 



(46) 



where e(p) is the polarization vector of neutral gauge boson. From Eq. (1461) . we get the 
counter terms at electroweak scale in on-shell scheme: 

,2 



m 



6Z° s (mJ=Al + J ^A* = A* + xA* 



5m 2 ' os (m ) = A Z U A 2 - m 2 5Z os 

z \ z/ U z z 



(47) 



Accordingly, the effective Lagrangian from the counter term diagrams is written as 
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3x, 



(1 

(a; + m;)] [(rf - Q f sl)\o; + o;) + qy^ot + 3 +) 



-QA( T f-QA) A l 



X* x c 



r (l + lnx l )l(Of- + 0: 



G> e x w T 2 (l + e) yJfcz^ 

v w w v ; >■ 



6v^' 

x 2 ^ 3 ' ; x 

2 



r 1 -^Ftv ~l~ x i 



2 

3x„ 



2 V3 V x 2 

z z 

( T f - QX) V; + o,-) + qv (0+ + o; 

2 2^> 11 (x F , x 



J 1/2 



ear 



; , ■ "Fp j 2 

+ — lnx z 



.I'- 



ll 



3x 2 



[(T/-gX)>7 + ^) + ^<W + ^ H 



1 ^ Fiy X l 



^ " X^ ( 2 +lnX ')J( 6 +0 6 

z z 

-6Q,(£ L £ fl + ^^)s 2 fT z -g fS 2N )(x F x F f' 2 



ex z 



61,1 \ X F a 5 X F n 



1 + In x, 



x' 



x- 



W + o:) + 



(48) 
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As a result of the preparation mentioned above, we use the contributions from the counter 
term diagram in FigJU(b) to cancel the ultraviolet divergence in Eq. (1421) . After matching 
the amplitude from FigH^a, b) to that from the diagram of effective theory (Fig^a)) and 
analyzing the corrections from light freedoms properly, we get the theoretical prediction on 
the lepton MDMs as 



%a+b) 



G F a e Q f m 2 x v 



R 1 2 



(\c L i 2 + ie 



.(49) 



24 v / 2tt 3 s 2 c 4 

w w 

In fact, the four-fermion interaction does not induce the corrections to lepton MDMs and 

EDMs in Figj5](a). Using the asymptotic expression of T 5 (x, z, u) at the limit z, u ^> x: 

1 

~9x 



T 5 (x,z,u) 



z + u . (z-u) 2 

-7T 2 [Z, U) H — 7T 3 [Z, U) + ■ ■ 



(50) 



Qx 18x 

one finally obtains the leading corrections contained in the asymptotic form of Eq. ([4*9!) under 
the assumption m F = m Fa = m F ^ ^> m w : 

ar (a+b) ^0(ml/m 2 F ) + -- - , (51) 

where ellipsis represents those relatively unimportant corrections. 
In a similar way, we can verify the identity 

kfA^(p 1 k) = 0. (52) 

A>° +d ) denotes the sum of amplitude for the diagrams in FigJH(c,d). After the steps adopted 
above, the corrections from FiglHc,d) to lepton MDMs and EDMs are formulated as 



l l(c+d) 



QpG F am 2 
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(53) 
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As a closed heavy fermion loop is inserted into neutral or charged gauge boson propagator, 
the counter terms to self energy diagrams of Z or W gauge bosons induce the renormalization 
for Weinberg angle 

5 s c 2 r5m 2 ' os 5m 2 ' os - - 2 



s,„ 2s s 



m? m? 



'-Ap F + ---, (54) 



2s 



where the dots indicate again nonleading contributions. Furthermore, Ap F denotes the 1- 
loop corrections from the heavy fermions to the p parameter, which appears in the ratio 
of weak neutral to charged current amplitudes and comparisons of the Z and masses. 
Using the above equation, one can express the corresponding counter terms for the vertices 
Zll as 



isa 



ie 



6s. 



+ sjs w u + . (55) 



Inserting the counter terms into one loop standard model diagrams, we finally obtain 

e 2 m 2 , „ o . „ a\Ss„ 
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V / www 
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G * m > 1 1 - 2s~ ■ ■ 2-s_' )\p, ■ .... m 



12 y/2- 



2 Q 2 



TT^S 



In principle, we should consider the corrections from the counter terms for W + uil: 

i5C w+Vl i = — -=— -5s w u- . (57) 

v w 

However, the corrections can be absorbed in the one-loop result as we parametrize the final 
result in terms of G F determined from the muon's lifetime. 

The contributions from the corresponding diagrams contain the additional suppressed 
factor m^/A^ w when both of virtual neutral gauge bosons in FigJH are replaced with the 
neutral Goldstone G°. Nevertheless, we should consider the corrections from those two loop 
diagrams in which one of virtual neutral gauge bosons is replaced with the neutral Goldstone 
G° since it represents the longitudinal component of charged gauge boson in nonlinear 
gauge. For many electroweak theories contain the neutral CP-even and CP-odd Higgs, we 
also generalize the result directly to the diagrams in which a closed heavy loop is attached 
on the virtual neutral gauge boson Z and neutral scalars h°, A . 

25 



(a) 



(b) 



FIG. 6: The two- loop diagrams in which a closed heavy fermion loop is attached to the virtual 
Z, G° or h° bosons. In concrete calculation, the contributions from those mirror diagrams should 
be included also. 



D. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual Z, G° (h°, A ) b osons 



Similarly, the renormalizable interaction among the electroweak neutral Gold- 
stone/Higgs G° (h°, A ) and the heavy fermions F a ^ can be expressed in a more universal 
form as 
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S°FF 
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)F p + h°F a (H s ag ^ + H £ ct lu J , 
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+ h.c. 



(58) 



where the concrete expressions of 7~l £ and T~t® depend on the models employed in 
our calculation, and the conservation of electric charge requires Qp = Q a - Generally, the 
couplings among the neutral Goldstone/Higgs and leptons are written as 



em, 



- iG°h 5 l + B £ h°ll - iB o A l~? 5 l 



(59) 



s °" 2m s 

W Vv 

and the parameters Bg, Bq depend on the concrete models adopted in our analysis. In full 
theory, the couplings in Eq. (|58l) induce the corrections to lepton MDMs and EDMs from the 
diagrams in Figj6j Accordingly, the corresponding diagram of effective theory is presented 
in Figj5]^b). In principle, we should also consider the corrections from those two loop dia- 
grams in which one of virtual Z gauge bosons is replaced with the neutral Goldstone/Higgs 
G° (h°, A ) in FigJU^a). However, the sum of amplitude from those two loop diagrams does 



not contribute to the lepton MDMs and EDMs when one of virtua 
FigHta) is replaced with the neutral Goldstone/Higgs G° (h°, A ) 



neutral gauge bosons in 



lit 
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After the steps taken in WW sector, we formulate the corresponding corrections to the 
lepton MDMs and EDMs from the diagrams in Figj6] as 

afk ° = - G tow!? {T ? ~ 2 ^/ S w)(^^) 1/2 [2(2 + \nx F0 ) Qo>1 (x„x h ) 

w w 



i T ? ~ 2Q f sl)(x w x F ) 1/2 [2{hix Fa - \nx Ffj )g 01 (x z ,x h ) 



128tt 3 s 2 c 2 

w w 

—F s (x x ,x h ,x Fa ,x Fg ) - F 6 (x z ,x h ,x Fg ,x Fa ) ^{yi^J,^ 

F e ; 



R 



G CX TYl^B) 



+F 5 (x. , x A , x Fa ,x Fg ) + F 6 (x, , x A , x F0 ,x F j] K(H° a ti L aB ~ 



df 



ZA o Gpa^mfio z 



(r; - 2Q f si ) (x w x F/j ) l ' 2 [2 (2 + In x F(j ) g 0A (x z , x A ) 



128vr 3 s 2 c 2 

w w 

+F 5 (x z , x A , x Fa ,x p + Kli, 



R 



ZG 



Gpa ^-{Tf - 2Q f sl)(x v x F0 )^[ - ^(\nx Fa - Ins, 



647T 3 S 2 c 2 

WW 

+F 5 (x. , * z , x Fq , x F0 ) + F 6 (x. , * z , x F/j ,x F j]U (f££ ~ CC ) 



(60) 



128tHs 2 c 2 

w v, 

+F 5 (x z ,x z ,x F ,x p )]$f(g"Z L B + gXi£R 

The expressions of form factors Fi(x,y,z,u) (i = 5, 6) can be found in appendix. In the 
heavy mass limit m F = vn Fci = m F0 3> m k , m A , m w , we have 
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d: 
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E. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual 7 and Z bosons 

When a closed fermion loop is attached to the virtual 7 and Z gauge bosons (FigJTj), 
the corresponding diagrams of effective theory are presented in FigJHJ Taking the steps 
above, one can get the tedious correction to the effective Lagrangian. If we ignore the terms 
which are proportional to the suppression factor 1 — 4s 2 , the correction from this sector to 
the lepton MDMs from this sector is drastically simplified as 



FIG. 7: The two- loop diagrams in which a closed heavy fermion loop is attached to the virtual 
7 and Z bosons. In concrete calculation, the contributions from those mirror diagrams should be 
included also. 



and the correction to the lepton EDMs is zero. In the limit m F = m Fa ^> m z , we approximate 
the correction to the lepton MDMs from this sector as 





(62) 




(63) 
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0") 



FIG. 8: The diagrams of effective theory corresponds to those diagrams in FigJTl Where the 
diagram (a) corresponds to Figf7|a) in full theory, and the diagram (b) corresponds to FiglT^b) in 
full theory respectively. 



F. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual 7, G° (h°, A ) bosons 

As a closed fermion loop is attached to the virtual neutral Higgs and photon fields, a 
real photon can be emitted from either the virtual lepton or the virtual charginos in the self 
energy diagram. When a real photon is emitted from the virtual charginos, the corresponding 
"triangle" diagrams belong to the typical two- loop Bar-Zee-type diagrams jig . Within the 
framework of CP violating MSSM, the contributions from two-loop Bar-Zee-type diagrams 
to the EDMs of those light fermions are discussed extensively in literature [7]. When a 
real photon is attached to the internal standard fermion, the correction from corresponding 
triangle diagram to the effective Lagrangian is zero because of the Furry theorem, this point 
is also verified through a strict analysis. Replacing the virtual neutral gauge boson with 
photon in FigEJ one obtains the relevant diagrams in full theory. Meanwhile, the diagram 
of effective theory is same as that presented in FigJSJ^a). After the steps adopted above, the 
corresponding corrections to lepton MDMs and EDMs from this sector are expressed as 

a l = ^3 1 ^( / ^f a )( X w X F a ) ^ ^11 ( X h J X F a 1 X F a ) 1 



8tt 3 



df = -^^£ 9 («fJ(x.x,jVV 12 (x„ 



d] A ° = - ^^^ ml)(^J 1/2 Tn( XA ,x r ,, Xr J . (64) 
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Note here that the corrections from this sector to the MDM of lepton depend on real parts 
of the effective couplings 7i £ and 7i° , and the corrections from this sector to the EDM 

1 ^ eta aa 1 

of lepton depend on imaginary parts of the effective couplings Ti.^ a and TC . In the limit 
m F = m Fa ^> m h , m A , the above expressions are simplified as 
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(65) 



F "A 

Similarly, the corrections to the lepton MDMs and EDMs from the jG sector are: 

,2 



jG 
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(66) 



The corrections from this sector to the MDM of lepton are proportional to real parts of 
the effective couplings , and the corrections from this sector to the EDM of lepton are 
proportional to imaginary parts of the effective couplings ', separately. In the limit 
m„ = m„ 3>m,we have 

G F oi e Q f Q a mfrn -< 2 



,-yG 
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167T 3 ?Ti I 
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1 + ln— F - 

m z 



(67) 



which are suppressed by the masses of heavy fermions. 



G. The corrections from the diagrams where a closed heavy fermion loop is at- 
tached to the virtual 7 bosons 

At the two loop level, there are QED diagrams involving a photon vacuumpolarization 
subdiagrams. For leptons or quarks, these contributions are of course known (3]. If heavy 
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fermions contribute, the two loop QED contributions are modified by the photon vacuum 
polarization (FigEJ). Adopting the same steps in WW sector, we formulate the corrections 
from Figj9]to lepton MDM and EDM respectively as: 



,77 



dp 



457T 3 







(68) 



which coincide with the well known results in Ref. 24|. For masses m Fa > lOOGeV, these 
corrections from this sector to the muon MDM a„ are below 10 -13 and hence negligible. 




W 7 

FIG. 9: The two-loop diagrams in which a closed heavy fermion loop is attached to the virtual 7 
bosons, where subdiagram (a) represents the two loop diagram in full theory, (b) represents the 
corresponding counter diagram, and (c) is the corresponding diagram of effective theory. 



III. THE LEADING TERMS IN TWO LOOP CORRECTIONS TO a t IN CON- 
CRETE ELECTROWEAK THEORIES 

In this section, we will present the leading terms from two loop corrections to a\ within 
some electroweak models discussed extensively in literature. 
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A. the SM 



Within the framework of SM, only top quark belongs to the fermion which mass is 
heavier than that of weak gauge bosons. The couplings in Eq.® and Eq. (j4"0l) are respectively 
expressed as 

>>tt 3 w 3 w 

C = . <£ = . (69) 

Using the unitary property of Cabibbo-Kabayashi-Maskawa (CKM) matrix V and m z = 
m w /c w , one formulates the leading corrections from top quark to lepton MDMs as 



s}l <7,.n. //r f 104 2 16 2 



m 



2 



SM = re , 3 _ ~- 2 _ ~ a ln 

2L Sv^vr^ 2 1 9 w 9 w m'j 

v w V h * 

-i^&t 1 - 2 ^ 2 ^- (70) 



In order to produce the terms oc m 2 , we take the limit s 2 — > 1/4 used in Ref4|, and 
approximate the last term as — 5G F m 2 Ap SM /24\/2ir 2 . Using the leading contributions 
Ap„.. = 3e 2 m 2 /4(47r) 2 s 2 m 2 + ••• in the limit m. 3> m , we recover the terms cx m 2 in 
Ref|4( perfectly. 

B. the Extension of SM with the Fourth Generation 

Besides top quark, the extension of SM with the fourth generation also includes addi- 
tional three heavy fermions: t', b', r 1 . Correspondingly, the couplings in Eq.(jHJ) and Eq. (|40p 
are separately written as 
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(71) 

- T , ■ 
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Assuming m , = m 



m 



m F ^> m w and applying the unitary property of 4 x 4 CKM 



matrix V, we finally obtain 
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Here, the 1-loop corrections to p-parameter from the heavy fermions of 4th generation can 
be written as 
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4G 16tts 2 m 2 

WW w 

In this model, the dominant contributions from Higgs sector are originated from the two- 
loop 7/1 diagrams. Under the assumption m , = m , = m , = m p , those contributions are 
zero since the anomalous cancelation. There is no correction from those heavy fermions to 
lepton EDMs also in the extension of SM with the fourth generation. 



C. the minimal supersymmetric extension of SM 



In this extension of SM, the additiona 
X12 an d four neutralinos (i = 1, ■ ■ • , 4) 
and heavy fermions are given as 

i L ± ± =2^cos2£ w + (f/t Uf/ J i ^ 

e ± . = 2^cos2^ w + (Ul) al (U R ) ig , (a,P = 1,2) 



possible heavy fermions include two charginos 



171 ] . The couplings among weak gauge bosons 







e o =^3^, (a,/? = l,--.,4), 
C L ^^(C/J^-^JVl^, 

C* ± =Ma(^ + -U»( l7 IU> (a = l,...,4; = 1, 2), (71) 

Xq-X^ v ^ 



Here, the 4x4 matrix M denotes the mixing matrix of neutralinos, two 2x2 matrices U L , U R 
denote the left- and right-handed mixing matrices of charginos, respectively. In the limit of 
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heavy masses, the mass spectra of charginos and neutralinos are respectively approached as 



m , ~ 

X 1 



m n ~ 



diag(|m 2 |, |// H |) , 
diag(|mi|, \m 2 \, \fi H \, \fi B \) 



(75) 



Here, p H represents the p parameter in superpotential, and m 2 , m\ denote the soft breaking 
masses of SU(2) x U(l) gauginos,respectively 

Applying Eq.( 1T4l) and Eq.(175l). we get the two loop corrections to lepton MDMs in the 
heavy mass limit \m\\ — |m 2 | 



\p H \ = m F > m w as 
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(76) 



the 1-loop corrections to p-parameter from the heavy supersymmetric fermions can be writ- 



ten as 
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(77) 



The abbreviation symbols eg = cos/5, sp = sin/5 with tan/3 = v 2 jv\ denoting the ratio 
between the absolute values of two vacuum expectations: V\, v 2 - As 1 < tan/? < 60 and 
m F ~ ITeV, the contributions from this sector to muon MDMs is well below 10~ n which 
can be ignored safely. 

At large tan/3, the dominant two-loop supersymmetric corrections to lepton anomalous 
dipole moments are originated from those Bar-Zee type diagrams which are analyzed ex- 
tensively. To obtain the corrections from those sectors, we formulate the relevant couplings 

as 
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(78) 



with Z z is the mixing matrix of two CP even Higgs. Assuming |/i H | = \m-i\ = \m\\ = m 
and 9i = 02 = Op, we expand the effective couplings in Eqj39]and Eql65]in powers of m w /m 
and get 
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where 0i i2 = arg(m 12 ), M = arg(/i H ) are the corresponding CP violating phases. Applying 
the above equations and Be = Bo = Be = tan /3, we find 
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Here, the form factors A, B, C, Z) depend on the masses of higgs and the mixing matrix 
of neutral CP-even higgs. 

The two loop corrections to lepton MDMs and EDMs are proportional to l/m 2 F in large 
tan/3 limit, which are consistent with the result presented in Ref.js| qualitatively. Using 



HME approximation and projection operator method, Ref.[5j approaches the corrections 
to the muon MDM from two loop diagrams in which a closed chargino/neutralino loop is 
inserted into those two Higgs doublet one loop diagrams as 

■jp-fls]) = ii x ir ^)( l -^y stgnM , (si, 

ou ivi slJSY 

under the assumption /i H = m 2 = (3c 2 /5s 2 )mi = m A = M SUSY and large tan/3 limit, where 
m A denotes the mass of neutral CP-odd Higgs. 
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M SUSY /GeV 

FIG. 10: The supersymmetric corrections to the muon MDM and varies with the super symmetric 
scale M SUSY when /i H = 771,2 = (^ c t/^ s t) m i = m A = -^susy an d tan/? = 5, 50. Here the solid 
line stands for the two-loop contributions from neutralino/chargino sector with tan/3 = 50, the 
dash line stands for the results of Eq. (l8ip with tan = 50; the dot line stands for the two- loop 
contributions from neutralino/chargino sector with tan/3 = 5, the dash-dot line stands for the 
results of Eq. (|8ip with tan = 5. 



To compare our result with that presented in Ref. [5J numerically, we take the same as- 



sumption /x H = m 2 = (3c 2 /5s 2 )mi 



m, 



SUSY 



on the parameter space of supersymme- 



try. In addition, the existence of a CP-even SM like Higgs with mass above 115 GeV sets a 
strong constraint on the parameter space of the employed model. To address this problem, we 
include all loop corrected effects in the Higgs potential 25] , and choose the Yukawa couplings 



of the 3rd generation sfermions as A t = A b = A T = M SUSY . Considering those points above, 
we present the numerical results in FigflQl For large tan /3 case, our numerical results agree 
with the approximation presented in Eq. (1811) very well. It implies that the equation Eq.( l8Tl) 
fits the exact result perfectly for large tan/3 and /i H = 777.2 = (3c 2 /5s 2 )mi = m A = M SVSY . 
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D. the littlest Higgs with T-parity 



In the framework of the littlest Higgs with T-parity, all the SM particles are even, as 
well as the corresponding mirror fields are odd under the discrete T-transformation [19]. In 
order to avoid dangerous contributions to the Higgs mass from one-loop quadratic diver- 
gences, we introduce additionally one T-even top quark T + together with its mirror partner, 
the T-odd top quark T_ besides the SM fermions /* and their mirror partners f l H . The cou- 
plings among the SM gauge bosons and heavy fermions in Eq. ( HOI and Eq. (jSJ) are respectively 
given as [20] 
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(82) 



Here, / is the breaking scale of a large SU(5)/SO(5) symmetry, and 



A 2 + A^ 



^ = 7^71 , (83) 



with A12 represent the Yukawa couplings of top quark sector. Additionally, the relations 
among the masses of heavy fermions are presented as 



T + v 
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m ^ H = m -* H ( 1 -^)' (^ = 1 ' 2 ' 3 )- ( 84 ) 

Applying the equations above, we give the leading corrections from heavy fermions to lepton 
MDMs in the limit f — m . = m . = m„ ^> m as 
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the 1-loop corrections to p-parameter from heavy fermions can be written as 2l|] 
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E. the universal extra dimension 
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If all particles of the SM are zero modes of corresponding 5-dimension bulk fields 22] , 
the KK excitations of fermion acquire the masses 



m /iw =^ +-> = !, 2, ■••), (87) 

where m. denotes the mass of corresponding SM field, and R is the compactification radius. 
To fit the present experimental data, we choose 1/R> 200 GeV. Furthermore, we formulate 
the couplings among the zero modes of weak gauge bosons and the KK excitations of fermions 
as [221 
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(88) 



with 
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(89) 



Using Eq. (!H71) . Eq.t lHBl and Eq. (l89l) . we formulate the leading contributions from the KK 
excitations of fermions as 



IV. CONCLUSION 

In this work, we have investigated the electroweak corrections to the lepton MDMs and 
EDMs from some two loop diagrams in which a closed heavy fermion loop is inserted into 
those two Higgs doublet diagrams. Adopting on-shell scheme, we subtract the ultraviolet 
divergence caused by the subdiagrams and get the theoretical predictions on lepton MDMs. 
As the masses of virtual fermions in inner loop are much heavier than the electroweak scale, 
we verify the final results satisfying the decoupling theorem explicitly if the interactions 
among Higgs and heavy fermions do not contain the nondecoupling couplings. Our results 
are universal for all extensions of the SM where the interactions among the electroweak 
gauge bosons and heavy fermions are renormalizable. As application of our analysis, we 
present the leading corrections to lepton MDMs in some popular extensions of the SM, such 
as the fourth generation, supersymmetry, universal extra dimension, and the littlest higgs 
with T-parity. 
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APPENDIX A: THE FORM FACTORS 
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APPENDIX B: THE FUNCTIONS 



The definition of *&(x, y,z) is written as: 

• A 2 > 0, + y/z < y/x: 

*(*, y,z) — 2 In + ln f-^ + ^~ A ) _ l n 1 l n £ 

v 2x ' ^ 2x ' 

-2L fa ( ' + »-'-y 2L fa (i^±i^) + f (Bl) 
where L i2 (x) is the spence function; 

• A 2 > 0, + ^ < ^y: 

*(a;,2/,«) = Eq.flEU)(z <->!/) ; (B2) 

• A 2 > 0, v/x + < y/z: 

y, z) = Eq.flEKa; ^ z ) ! ( B3 ) 
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• A 2 < 0: 
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where Cl 2 (x) denotes the Clausen function. 
The expressions of <p (x,y), ipi(x,y), ip 2 (x,y) and <^ 3 (x, y) are given as 
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The functions adopted in the text are written as 
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